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An overview is given of the study of final state interactions in hadronically decaying W pairs 
produced in e^e^-collisions as it is performed by the four LEP experiments. Bose-Einstein 
correlations are investigated by comparing like- with unlike-signed pairs of pions and/or using 
the mixed event analysis technique. Colour reconnection is examined with a method that 
compares the particle flow distributions in inter-jet regions. 



1 Introduction 

The second phase of the LEP programme was devoted to run at centre-of-mass energies above 
the W pair threshold. When such a pair is produced, the W bosons will immediately decay due 
to their very short lifetime (r^y ~ 0.1 fm/c). In 46% of the cases this will give rise to a final 
state with four quarks, the so-called hadronic channel. Considering the full integrated luminosity 
delivered by the LEP machine, each experiment has collected roughly 5000 such events. The 
produced quarks will hadronise after typical scales of about 1 fm, resulting in jets of particles 
measured with the detector. Due to the large space-time overlap, correlations between the 
particles coming from different W's might occur. The possible effects studied at LEP are Bose- 
Einstein correlations (BEC) between identical bosons in the jets, and colour reconnection (CR) 
between the quarks themselves. 

When an exchange of particles or energy and momenta between the two W bosons is in- 
adequately described by the Monte Carlo programs, this will cause a bias in the measurement 
of the W mass. The statistical uncertainty on the W mass is 30 GeV/c^ when combining the 
results from all the LEP experiments using the hadronic channel onlyS. Various models exist to 
describe the Final State Interactions (FSI) and a number of them predict biases several times 
larger than this statistical error. Besides the measurement on itself, these results are therefore 
also important for the determination of the W mass. 



Table 1: The energy ranges and corresponding integrated luminosities used by the four LEP experiments for the 
studies of Bose- Einstein correlations and colour reconnection. 



Experiment 


BEG 

^ [GeV] C [pb-i] 


CR 

^ [GeV] C [pb-i] 


ALEPH 
DELPHI 
L3 
OPAL 


183 - 209 683 
189 - 209 531 
189 - 209 627 
172 - 189 250 


189 - 209 628 
183 - 209 601 
189 - 209 627 
189 183 



2 Bose-Einstein Correlations 

Since the bosons produced in the hadronisation process obey Bose-Einstein statistics, an en- 
hancement in the number of identical bosons created close in phase space is expected due to 
amplitude symmetrisation. This effect has been observed in e.g. Z ^ qq at LEPl, thus be- 
tween particles coming from the same boson. Now that boson pair production is possible at 
LEP2 energies, the question which arises is whether BEG also exist between particles coming 
from different W's (inter-W BEGL The energy ranges and luminosities used by the four LEP 
experiments for this measurement B can be found in table |l|. 

BEG are typically studied by looking at the two-particle density in momentum space 

where Upairs is the number of pairs in N^,^ selected events and Q is the 4-momentum difference 
\/ —{pi — P2)^^ with pi and p2 the momenta of the two particles. The presence of BEG will give 
rise to an enhancement of like-signed particle pairs at low Q. The two-particle density is often 
parametrised as p ~ (1 + Aexp(— iJ^Q^)) and interpreted as a Gaussian source of radius R and 
strength A. 

The BE32 algorithm from the LUBOEI model i is used to generate Monte Garlo samples 
containing BEG inside and between the W bosons. This model just reshuffles the momenta of 
identical bosons and is tuned to the hadronic Z decay data to describe BEG in these events. 



2.1 Unlike-sign Pair Analysis 

In this method, the ratio between the number of like-signed pairs and unlike-signed pairs ob- 
tained from data is compared with a Monte Garlo sample generated without inter-W BEG and 
studied clS cl function of Q. 

OPAL extracts these values for three event classes: hadronic and semi-leptonic WW decays 
and non-radiative Z* — > qq events. The double ratios are then deconvoluted to obtain the 
contributions of pions coming from different W's, same W and hadronic Z decays. The result is 
shown in figure ||. A fit to these distributions results in a strength parameter A = 0.69ib0.12ib0.06 
for pairs coming from the same W, again proving the existence of BEG inside the W. The fit 
to the distribution for pions from different W's gives A = 0.05 it 0.67 it 0.35, being compatible 
with both inter-W and no inter-W BEG due to the small statistical sample used. 

The double ratio from the ALEPH data up to 202 GeV shown in the right plot of figure 
|l], is fitted with a similar function as described above and the integral over Q of this function 
is calculated. The same procedure is followed for a MG sample using the JETSET model with 
inter-W BEG. When compared to the data, this model is disfavoured at the level of 2.2a. 
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Figure 1: The double ratio distributions as a function of Q from OPAL (left), deconvoluted to different contribu- 
tions, and ALEPH (right), where the fit to the data is indicated as a solid line. 



2.2 Mixed Event Analysis 

Inter-W BEC can be accessed directly from data, without relying on a certain model, by com- 
paring the hadronic decays with events constructed from the hadronic part of two different 
semi-leptonic events This analysis has been preformed by ALEPH, DELPHI and L3. 

If the two W bosons in a hadronic event decay independently, the two-particle density can 
be written as 



P 



2/ 



(2) 



where the first term on the right hand side is taken from semi-leptonic events, and the second 
term is obtained by mixing the hadronic part of two different semi-leptonic events. Then the 
sensitive distributions are the difference or the ratio between left and right hand side: 
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When BEC between particles coming from different W bosons do not exist, Ap = and D = \. 
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Figure 2: The Ap distribution as a function of Q from ALEPH (left) and the integral J = 

different centre-of-mass energies from L3 (right). 
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Figure y shows the Ap distribution for the ALEPH data, which is in good agreement with 
the prediction from a MC containing only intra-W BEC. The same conclusions can be drawn 
when looking at these distributions from DELPHI and L3. Another way to look at the data is 
to calculate the integral J = /o^"""" Ap{Q)dQ. As can be seen in the right plot of figure |2|, the 
value obtained from a MC with inter-W BEC clearly does not agree with the data. 

Comparing the strength parameter A from 
the fit to the distribution of the ratio D from 
the DELPHI experiment (fig. the MC con- 
taining inter-W BEC is excluded on the level of 
3.2(7. L3 looks at the ratio between the D dis- 
tribution from data and this distribution from a 
MC with only intra-W BEC in order to eliminate 
some possible systematic errors. Comparing the 
strength parameter from the fit to this double 
ratio using the data with the one obtained when 
a MC with inter-W BEC is used, the latter one 
can be excluded by 4.7 standard deviations. 
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Figure 3: The D distribution as a function of Q from 
the DELPHI experiment 



3 Colour Reconnection 

In the string picture, the hadronisation particles come from the decay of the colour string 
stretched between the quarks coming from the same W. When the colour flow pattern is modified, 
referred to as colour reconnection, strings might span between two quarks from different W 
bosons and sebsequently decaying particles cannot be uniquely assigned to either W. 

Different phonomenological models exist to describe CR. In the model of Sjostrand-Khozei, 
as it is implemented in PYTHIA, CR might occur when strings overlap. In the type I model 
(SKI), the strings are associated to colour flux tubes and the reconnection probability, Preco, is 
related to the volume of the overlap, f{y/s), by the following relation: 

P.eco = 1 - e-^(^)•^^ (4) 

where kj is a free parameter that can be varied to obtain samples with different reconnection 
probabilities. Also the predictions from other models like SKII, ARIADNE, HERWIG and 
Rathsman are tested against the LEP data. 



3.1 Particle Flow Method 

This method compares the particle rates between jets from same and different W bosons, since 
CR will result in a depletion or enhancement in inter-jet regions. This approach, pioneered by 
Lsi, is much more sensitive to CR effects than previous methods and is now adopted by th£ 
other experiments. The data samples used by the different LEP experiments for this analysisu 
are indicated in table |^. 

In order to study the particle flow, the momenta of the particles are projected onto the plane 
spanned by the most energetic jet and the jet most likely coming from the same W. Since the 
W decay planes are not planar, the momentum of a particle i is further projected on the planes 
spanned by the 2 adjacent jets k and j. In order to be able to compare the different inter-jet 
regions, the angle (pi between the projected momentum and the most energetic jet, is rescaled 
by the angle (p^j between the nearest jets. 

DELPHI and L3 perform a topological selection with tight cuts on the angles of the jets in 
order to obtain a high purity sample with almost planar jets, but with a low selection efficiency. 
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ALEPH and OPAL use a selection based on there standard W mass selection, which has a larger 
efficiency but the jet topology is less trivial and the background contamination is higher. 

An example of the particle flow distribution 
is shown in figure |^. Due to the definition of the 
rescaled angle, the direction of the jets are point- 
ing to integer values. The two regions inside a 
W (indicated with A and C on the plot) are now 
combined, as well as the two regions between the 
W's (B and D). These two particle flow distribu- 
tions are then compared to each other by taking 
the ratio, i?, which can be seen in the left plot 
of flgure |5|. In order to quantify the CR effects 
of the different models, the particle flow inside 
and between the W's is integrated over (\>resc be- 
tween 0.2 and 0.8, the region most sensitive to 
CR. Then the ratio between the two integrals is 
called i?Ar and shown in figure ^. 




Figure 4: The particle flow distribution as a function 
of the rescaled angle from DELPHI. 
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Figure 5: The ratio R between the particle flows inside and between the W's (left) and the ratio Rn between the 
integrals over the particle flows as function of centre-of-mass energy (right). 



In order to compare the different experiments, the ratios between the i?7v value from data 
and a MC without and with CR (SKI with 100% reconnection probability) have been calculated! 
and are shown in table ||. A deviation from unity is a measure of the disagreement between data 
and the model used. The measurements of DELPHI and L3 show clearly a preference for no CR, 
while ALEPH prefers a model inbetween no and 100% CR. The OPAL experiment comes to the 
same conclusion when using the W mass selection. They also applied the topological selection 
as a cross check analysis and then found a preference for no CR. This is actually under study in 
all LEP experiments to see whether this is due to the type of selection, the background or other 
possible sources. 

The results can now be interpreted in terms of the free parameter ki in the SKI model, see 
equation (Q), by calculating a from the comparison of Rn between data and MC samples 
evaluated at different values of kj. L3 obtains a limit of kj < 1.55 at 68% C.L., while ALEPH 
prefers a value of (kj) = 3.55 (a reconnection probability at y/s = 189 GeV of about 50% and 
80% respectively). 



Table 2: The ratios between the Rm variables from data and a MC sample, one without CR and one with CR 
using the SKI model with 100% reconnection probability. The values between brackets indicate the deviation 

from unity. 





{Rr) Data/no CR 


{Rn) Data/CR 


ALEPH 
DELPHI 
L3 
OPAL 


0.961 ±0.012 ±0.007 {-2.8a) 
1.009 ±0.030 ±0.019 (±0.3ct) 
0.990 ± 0.025 ± 0.023 {-0.3a) 
0.906 ±0.033 ±0.011 (-2.7a) 
0.996 ±0.051 ±0.011 {-0.1a) 


1.041 ±0.013 ±0.008 (±2.7cr) 
1.110 ±0.033 ±0.029 (±2.5cr) 
1.194 ±0.039 ±0.029 (±4.7o-) 
1.050 ±0.038 ±0.013 (±1.2cr) 
1.193 ±0.061 ±0.014 (±3.1cr) 



4 Conclusions 



A large part of the data samples collected by the four LEP experiments has been analysed 
to study possible final state interactions in hadronically decaying W pair events. None of the 
experiments observe any indication for the existence of Bose-Einstein correlation between par- 
ticles coming from different W bosons. The LUBOEI model used to simulate such an effect is 
strongly disfavoured. The results from the different colour reconnection analyses are inconsis- 
tent. DELPHI and L3 prefer no CR, while ALEPH prefers a large value for the reconnection 
probability. OPAL gets different results depending on the selection used. All experiments are 
presently investigating the possible causes. When the foreseen combination of the results from 
the four collaborations will be performed, the contribution to the systematic error on the W 
mass from the effects described here will be reduced significantly. 
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